arXiv: 1507.03008vl [astro-ph.SR] 10Jul2015 


Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 13 July 2015 (MN style file v2.2) 


Infrared Variability from Circumbinary Disc Temperature 
Modulations 

Eva H.L. Bodman, & Alice Quillen 

Department of Physics and Astronomy, University of Rochester, Rochester, NY 14627, USA; 


ABSTRACT 

The temperature of a circumbinary disc edge should undulate due to variations in 
illumination as a function of binary orbital phase. We explore circumbinary disc tem¬ 
perature variations as a source of broad-band infrared light curve variability. Approx¬ 
imating the wall of a circumbinary disc edge as a wide optically thick cylinder with 
surface temperature dependent on its illumination, we find that a binary comprised 
of 1 Mq and 0.5 Mq pre-main sequence stars in a ~15.5 day period, would exhibit 
the largest amplitude variations of ~9% at 3.77 and 4.68 ^m as seen by a distant 
observer. The amplitude of variations and shape of the light curve is sensitive to the 
luminosity and mass ratios of the stars in the binary, the radius of the circumbinary 
disc clearing, the binary separation, and the orbital inclination. The light curve varia¬ 
tions are smooth and very red with a non-sinusoidal shape for most of the parameter 
space explored. Possible morphologies include a single peak with a flat region, two 
peaks of different heights or a single dip. 
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1 INTRODUCTION 

Variability is ubiquitous among young stellar objects 
(YSOs) and occurs on a variety of time-scales and over 
a wide range of wavelengths, from radio through X-rays. 
Near and mid-infrared monitoring campaigns probe rela¬ 
tively cool material that is dim at shorter wavelengths, 
and are sensitive to distribution, temperature, and opacity 
changes in circumstellar gas and dust. Periodic variations 
in photometric measurements can be associated with rota¬ 
tional modulation of cool star spots or hot accretion spots, 
stellar pulsation, transits, eclipses and intervening circnm- 
stellar material (e.g., Fernandez & Eiroa 1996; Carpenter 
et al. 2001; Winn et al. 2004; van Boekel et al. 2010; Morales- 
Calderon et al. 2011; Wolk et al. 2013; Cody et al. 2014; 
Parks et al. 2014). 

Most young stars are born in binary star systems 
(Duchene & Kraus 2013) and in these, additional processes 
contribute to periodic light curve modulations including 
eclipses, ellipsoidal tidal deformations of the stars, reflec¬ 
tion of light of each star from its companion and at low 
amplitude, relativistic beaming (Zucker et al. 2007; Tal-Or 
et al. 2015). A binary star system can be embedded in a 
circumbinary disc (e.g., Jensen et al. 2007; Ireland & Kraus 
2008; Nagel et al. 2012; Biller et al. 2012; Takakuwa et al. 
2013; Gillen et al. 2014). Periodic broad-band photomet¬ 
ric variations could constrain the location or shape of the 
circumbinary disc edge (Nagel et al. 2010, 2012), infer the 


presence of a low mass companion (Demidova et al. 2013; 
Maiz Apellaniz et al. 2015) or discover a circumbinary or 
circumsecondary disc (Rattenbury et al. 2015). As studies 
increase light curve quality, the number of data points, and 
the wavelength coverage, it may become possible to differ¬ 
entiate between different variability mechanisms using light 
curve morphology and colour. Circumbinary discs are of par¬ 
ticular interest as recently discovered circumbinary planets 
imply that circumbinary discs host planet formation (Mar¬ 
tin & Triaud 2015). 

A possible source of IR variability is temperature os¬ 
cillations in a circumbinary disc due to uneven irradiation 
from the binary star system (see Nagel et al. 2010). The 
circumbinary disc centre is expected to be at the binary’s 
centre of mass for circular orbits. The disc can contain an 
inner clearing or cavity caused by the binary (Artymowicz 
& Lubow 1994; Pichardo et al. 2005, 2008). Since neither 
star in a binary is located at the centre of the disc, the dis¬ 
tance between one of the stars and a point in the inner edge 
or wall of a circumbinary disc depends on the binary’s or¬ 
bital phase. The oscillations in distance cause a modulation 
in the absorbed stellar flux locally in the disc which leads 
to a local disc wall temperature dependent on binary or¬ 
bital phase and azimuthal angle (see Figure 1). Nagel et al. 
(2010) studied the effects of a binary star on the tempera¬ 
ture of the inner wall of an optically thick circumbinary disc 
and the impact on the spectral energy distribution (SED) 
for the binary transition disc system CoKu Tau/4. The au- 
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thors predict variability in the SED with the orbital period, 
however, because of the long orbital period (about 20 years) 
for CoKu Tau/4’s binary they did not predict the shape of 
the light curve. Shadmehri & Khajenabi (2015) and Clan¬ 
ton (2013) investigated the effects of a binary star on the 
ice line in a circumbinary disc and focused on wider bina¬ 
ries. In both cases, temperatnre variations are due to the 
oscillations in distance from each star to the disc wall. 

In this paper, we explore how temperature variations 
in a circumbinary disc would introduce variability in a light 
curve. Our goal here is to estimate the amplitude, shape 
and colour of variations in a light curve that are caused by 
modulation in the temperature of an optically thick wall of 
a circumbinary disc. In Section 2, we describe our model 
and how we compute infrared light curves. In Section 3, we 
discuss the amplitude of the magnitude variations predicted 
by our model and how the amplitude is effected by the stellar 
properties of the central star and the geometry of the disc. 
We describe the sensitivity of the light curves to the stellar 
luminosity and mass ratios, orbital separation, circumbinary 
disc wall radius and orientation or viewing inclination. A 
summary and discussion follows in Section 4. 


2 THE MODEL 

We first describe the geometry of our model and then we 
describe how we integrate the flux from the disc wall. Our 
simple model consists of a circumbinary disc wall and two 
stars in circular orbits (see Figure 1). We assume that the 
disc is not inclined with respect to the orbital plane. The 
disc wall is a cylinder characterised by its vertical height, 
H, measured from the orbital mid-plane and its radius, Rd 
measured from the centre of mass of the binary. We make 
no assumptions about the composition of the disc wall other 
than it is optically thick and is well described as a blackbody 
with temperature dependent on azimuthal angle, 0 (see Fig¬ 
ure 2), neglecting temperature variations as a function of 
height above the orbital plane. We are neglecting dust fea¬ 
tures and only using a blackbody spectrum to describe the 
disc wall emissions. The dust features are expected to change 
the size of the light curve variations only by a factor of a cou¬ 
ple (see Section 3.5). We do not expect an optically thin disc 
to exhibit light curve modulations from temperature varia¬ 
tions as the total emitted disc flux would be independent 
of viewing angle. We assume the stars are approximated 
with blackbodies and that the binary orbit and circumbi¬ 
nary disc wall are circular. Both stars irradiate the disc but 
since neither star is at the centre of the disc, the distance 
between star and different points along the wall varies so 
the absorbed flux and thus the local temperature varies as a 
function of azimuthal angle and orbital phase, (j) = nt with 
(j) = 0 with stars oriented on a line perpendicular to the line 
of sight and n the mean motion of the binary. At (() = 0, 
the primary star is located at 6^ = 0 and the secondary star 
at 0 = TT so the primary and secondary eclipses, if observ¬ 
able, would occur when (j> = 7r/2 and 37r/2, respectively. As 
the stars travel along their orbits, the distance between each 
star and a point on the disc wall varies causing the observed 
flux from that region to oscillate. The half of the disc wall 
that contributes to the observed disc flux is the range 6 = 0 
to TT. 



Figure 1. Diagram of the temperature variation along a cir¬ 
cumbinary disc wall as it varies during the binary orbit. The disc 
wall is shown as a cylinder with colour showing its local tem¬ 
perature. Yellow is hotter than red which is hotter than brown. 
Orbital phase is denoted by </>. The disc wall is assumed to be 
optically thick so we have neglected radiation from its back side 
(shown in navy blue). The disc outside the wall is not drawn. At 
(p = ±7r/2, one star can shadow the other giving a cooler (brown) 
region on the disc wall. The primary, larger star (shown in yellow) 
is more luminous and the disc wall is hotter near this star. The 
observed flux at <)> = tt is the same as at <)> = 37r/2 but that at 
= 0 is not the same as that at 0 = tt. 

Each star is characterised by its effective temperature, 
radius, Ri, mass, Mi, luminosity, and semi-major 
axis, fli, with index i = 1,2 depending on the star. Their 
SEDs are given by blackbody curves. The separation be¬ 
tween the two stars we denote Usep- We neglect stellar 
eclipses, ellipsoidal deformation, reflection and relativistic 
beaming. We assume that the disc wall reaches thermal equi¬ 
librium on a time-scale much shorter than the orbital period. 
Following Nagel et al. (2010), the thermal time-scale is 

where k is the Boltzmann constant, mn is the mass of 
Hydrogen, a is the Stefan-Boltzmann constant, T is the 
temperature of the wall. The surface density of the wall is 
Efl = and the volumetric density is estimated 

hy p = "Ed/H. Note we are interested in the surface parallel 
to the disc wall and p is estimated with the surface density 
of the top of the disc. The surface density of a disc at about 
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Figure 2. Schematic top-down diagram of binary system with 
circumbinary disc illustrating the geometry of our model. The 
primary and secondary stars are labelled with mass Mi, M 2 , and 
have semi-major axes ai, a 2 and zero eccentricity. The circumbi¬ 
nary disc wall radius is Rd . The azimuthal angle 0 gives a location 
in the disc wall with 0 = 0 normal to the line of sight. The or¬ 
bital phase is described with (p = nt with 0 = 0 when the binary 
is perpendicular to the line of sight. Here n is the mean motion 
of the binary. Xi{9,(p) and X2{0,(p) refer to distances between 
a location in the disc wall and the primary and secondary stars, 
respectively. The line of sight is towards the top of the page. We 
define an angle 7 = 0 — 0 that helps us compute the distances 
Xi,X2. 


0.1 AU is Ed 500 g cm“^ (D’Alessio et al. 1998). Since the 
wall emission come from a thin radial region of the disc, we 
estimate the wall surface density by setting AR = H such 
that Ej? = Ed. The thermal time-scale is then 

The thermal time-scale is comparable to the orbital time- 
scale only for very close binaries. 

The visible part of the disc wall is the side facing the 
observer and the other side is assumed to be completely oc¬ 
culted by absorption of the outer regions of the disc (see Fig¬ 
ure 1). We model the disc wall as an optically thick cylinder 
with vertical walls and we neglect radiative transfer effects. 
All of the stellar radiation that reaches the wall is absorbed 
and reemitted at the surface. The disc wall is the warmest 
part of the disc and it receives a sizeable portion of the 
stellar irradiation. Consequently emission from the disc wall 
dominates the spectral energy distribution and we neglect 
emission from the rest of the disc. 

Assuming the disc wall has an albedo of zero, the ab¬ 
sorbed stellar flux at a given azimuthal angle 0 on the disc 
wall is 


^Wall(^, D = 


T* 


-b 


^nXU0,<l>) 


(3) 


when one star is not blocking light from the other star to 
the disc wall. The absorbed flux also depends on the orbital 



Figure 3. Diagram of the shadowed regions from the secondary 
eclipse on top and primary eclipse on bottom. The colours vari¬ 
ance along the disc wall indicates the local temperature in the 
same manner as in Figure 1. Yellow is hotter than red and navy 
is the back side where radiation is neglected. The brown colour 
shows the cooler temperatures from the stellar eclipses. The larger 
yellow star is the primary and the smaller red star is secondary. 
The lines mark the geometry used to calculate the width of the 
shadow which ignores the radius of the smaller star. The top dia¬ 
gram shows the shadow from the secondary eclipse (see equation 
5) and the bottom diagram shows the primary eclipse shadow 
(see equation 6 ). This approximation gives the shadowed region 
sharp edges in the centre of the smooth transitions shown in the 
diagram. The disc outside the wall is not shown. 


phase of the binary as the distances to the stars also depend 
on this angle. Here Xi{6, 0) is the distance between the i-th 
star and the point on the wall with azimuthal angle 0 in the 
orbital plane (see Figure 2) and during orbital phase 0. It is 
convenient to define an angle 

'y = 0-(f> (4) 

After choosing the semi-major axes, the distances Xi depend 
only on the angle 7 so henceforth we write Xi{'y). We use 
small angle approximations so that the flux is independent 
of height and the temperature variation as a function of 
height can be neglected. This is valid when H Rd- We 
also assume that the stellar flux is normal to the surface of 
the wall and this is valid when radii <C Xi. 

When one stars occults the other star, as seen from the 
disc, the disc temperature can be reduced (see Figure 3). 
Since the disc surrounds the binary star and is aligned with 
the binary orbit, there will be two regions along the disc wall 
where one stars shadows the other, as seen from a particular 
position on the disc wall. The shadows cause sharp dips in 
temperature and they are illustrated as brown regions in 
Figure 3. The primary’s shadow occurs at 7 = 0 while the 
secondary’s is at 7 = tt. We first discuss the angular width of 
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the primary shadowing the secondary and then the angnlar 
width of the secondary shadowing the primary. 

We first consider the sitnation where the primary is 
larger than the secondary. In this case, light from the 
secondary does not reach the disc wall during secondary 
eclipses. The shadow on the disc has an angular width Ay 
measured from the centre of the shadow where 

Itanyl ^ tan(A7) = ^ + a2 

-tXD Ojsep 

as illustrated by the top drawing in Figure 3. We denote 
the separation between the stars as Osep- The shadow where 
light from the primary is blocked by the secondary has a 
width of 

|tan(7r-7)| s; tan(A7) = (6) 

<^sep 

In this region the primary is only partially occulted by the 
smaller secondary and the shadow is smaller than the pri¬ 
mary’s shadow (see bottom of Figure 3). 

If the secondary is larger than the primary,^ then the 
primary partially occults the secondary where 

Itanyl sj tan(A7) = —— ( 7 ) 

-ttD ^sep 

and the primary is completely shadowed by the secondary 
with 

|tan(7r - y)] ig tan(A7) 

<3;sep 


In Figure 3, the width of the shadow is larger than the 
height of the disc but for stars with small radii or for high 
disc walls, that is not true. Here we neglect the height of 
the shadows, assuming that the shadow covers the full disc 
height for all star radii and disc heights. Since we are inte¬ 
grating over an area much larger than the shadow, the tran¬ 
sition or partially shadowed regions can be neglected. We 
assume that shadow region begins when the half of one star 
occults the other giving a sharp transition from no shadow 
to full shadow. 

For angle y satisfying equations 5 and 8, the disc wall 
only receives light from the primary or the secondary, re¬ 
spectively. In these shadows, the absorbed stellar flux is 


Twall.c 




where £ refers to whichever star is larger. For angles satisfy¬ 
ing equations 6 or 7, the disc wall absorbs a flux 


F = 

W-ll.P 47rXa(7)2 

Li f fRo + aiY] 

~^ 47 rXi(jy[ Rj{RD-asJ ) 


( 9 ) 


where s denotes the smaller star. Since the smaller star is 
much closer to the disk edge than the larger star, the smaller 
can completely occult the larger star when Re < RsX {Rd + 
ae)/{RD — Os)- In this case, the flux is the same as Twaii.c 
but with s replacing £. 


^ This does not occur for coevolving PMS stars studied here but 
may occur for more evolved systems. 


Assuming that the absorbed flux is reradiated ther¬ 
mally, the local mid-plane disc wall temperature outside of 
shadows (none of the conditions in equations 5 - 7 is satis- 
hed) is given by 


T(7) 


/ fwall (y) \ 

) 


5778 K 




( 10 ) 


where Fq = Lq/AtvRq is the surface flux of the Sun, 
is in solar luminosities and Xi is in solar radii and we have 
assumed a wall emissivity of 1. For regions of complete and 
partial shadow, Fwaii is replaced with Fwaii.c and Twaii.p, 
respectively. 

The azimuthal temperature profile has two peaks cor¬ 
responding to the regions along the disc wall closest to the 
primary at y = 0 and the secondary at y = tt. The relative 
sizes of the peaks depend on the primary to secondary lu¬ 
minosity ratio, L*,i/L*, 2 . The closer the ratio is to one, the 
closer to equal height the peaks are with the primary peak 
higher. For two identical stars the peaks at y = 0, tt would 
be the same height and temperature oscillations would oc¬ 
cur at twice the orbital period along with the light curve 
modulations. The secondary peak disappears at large lumi¬ 
nosity ratios as the primary dominates the stellar irradiation 
of the disc wall. The primary peak is lower for very large pri¬ 
mary to secondary mass ratios, M 1 /M 2 . When the primary 
is much more massive than the secondary, its orbit is small 
(from the centre of mass) so that its distance to the disc wall 
does not vary much. Sharp dips in temperature occur in the 
centre of each peak from each star shadowing the other and 
these occur at y = 0, tt. 

To model the observed light curve, the local disc flux 
is calculated using Planck’s function Bx{T{9,(j>), X) and 
the local disc wall temperature. We compute the observed 
monochromatic flux from the disc wall as a function of wave¬ 
length, binary phase cf>, and orbital inclination i by integrat¬ 
ing the local surface flux over the disc wall area 


mi)xDisc = ^ £ Bx{T{e,<P),X) 

xHeg{6,i) smisinOdd (11) 


where d is the distance to the disc from the observer and 
inclination i with an edge-on orbit having i = 90°. HeniO, i) 
is the effective disc height taking into account occultation of 
the wall by the nearer side of the disc as seen by the observer. 
Since the disc wall is optically thick, we assume the no ra¬ 
diation from disc wall reaches the observer in the regions 
obscured by the opposing side of the disc. The amount the 
nearer half of the disc obscures the wall depends not only on 
the inclination but also the azimuthal angle along the wall, 
6. The height of the projected hole in the disc is smaller near 
0 = 0, TT so more of the inner disc wall will be occulted. The 
disc is locally occulting itself if 

Rd cot isinO < H. (12) 


When true. Bed = Rd cot i sin 0 but otherwise the full disc 
wall height H is used. This self-occultation of the disc only 
gives a large effect at inclinations very close to 90° since 
we assume that the disc radius is much larger than the disc 
scale height. Lastly, we add the flux from each star to the 
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Table 1. Fiduciary Model Parameters 


Ml 

1 .OM0 

M2 

0.5Mq 

Ri 

1.935 Rq 

R2 

1.461 Rq 

Teff.l 

4417K 

Teff.Z 

3932K 

^^sep 

30 Rq Ri0.14AU 

i 

75 deg 

Rd/^ sep 

1.8 

H/Rd 

0.0960 


Fiducial model parameters for temperature profile and light 
curves in figure 4. Stellar mass, radius and surface temperature 
are consistent with 2 Myr old pre-main-sequences stellar 
evolution models by Siess et al. (2000). 


disc flux to calculate the total flux for the light curve, 

F{4>, t)A, total = i)^.Disc + Fx,! + Fa.2 . (13) 

Since the disc flux increases with disc temperature, there are 
peaks in the light curve at binary phases that give peaks in 
the temperature at 9 = -k 12. The light curve peak from the 
primary occurs at binary phase (j> = 7r/2. If the secondary 
temperature peak is large enough, there is also peak in the 
light curve peak at binary phase (f) = 37r/2. 

The shadows from primary and secondary typically 
have a small effect on the light curve. The shadows change 
the shape of the light curve slightly near = 7r/2 and 37r/2. 
The regions of lower temperature occur at the peaks in the 
temperature and decrease the observed disc wall flux the 
most when the temperature peaks are at 9 = -n 12. This de¬ 
creases the peaks in the light curve and increases dips that 
occur at (j> = 1^/2., 37 r/2. Occultation from the near disc wall 
does not change the shape but deceases the light curve peaks 
at very high inclinations. 

Three additional sources of occultation are neglected: 
the star occulting the disc wall, the disc occulting the star, 
and the stellar eclipses (one star eclipsing the other). The 
disc eclipsing the wall and the wall eclipsing the star occur 
approximately at inclinations larger than 80° but the exact 
value depends on the geometry of the system like the disc 
radius and the binary separation. Stellar eclipses typically 
occur at higher inclinations than the disc eclipses. The de¬ 
crease in observed flux from stellar eclipse or disc occulting 
the star is typically much larger than the variations from 
temperature oscillations since they are of order of the stel¬ 
lar flux, see Table 3. 


3 FIDUCIAL MODEL 

We begin our study of the light curve with a fiducial model 
and then we vary each model parameter from the fiducial 
model values individually. Our model is an approximate rep¬ 
resentation of a binary star with an optically thick circumbi¬ 
nary disc and such discs can be found around PMS stars. 
Young PMS stars typically have an optical thick inner region 
in their proto-planetary discs that extend outward from the 
star to about 10 AU (Williams & Cieza 2011). We choose a 
fiducial model consisting of two PMS stars of the same age 


and a circumbinary disc. We set the age to 2 Myrs which is 
well within the disc half-life of 5-6 Myrs (Bell et al. 2013). 
We use PMS stellar evolution models from Siess et al. (2000) 
to estimate the properties of a 1.0 Mq primary and 0.5 Mq 
secondary. We choose a binary separation of 30 Rq «0.14 
AU which sets the orbital period to about 15.5 days. The 
inclination is 75° which is low enough for neither star to 
eclipse the disc or the disc to eclipse either star. The val¬ 
ues for the stellar and orbital properties used in our fiducial 
model are listed in table 1. 

The disc radius and height can be constrained with dy¬ 
namical arguments. For a circumbinary disc, a hole in the 
centre is created through dynamical interactions with the 
binary. Artymowicz & Lubow (1994) found the minimum 
radius for the wall of a circumbinary disc through a hydro- 
dynamical approach that is proportional to the binary sep¬ 
aration and increases with eccentricity. For zero eccentricity 
orbits, the authors found Rd = l.Sttsep which we use for our 
fiducial model. With a separation of 30 Rq, the inner disc ra¬ 
dius is 54 Rq «0.25 AU. We estimated the height of the disc 
wall from hydrostatic equilibrium condition, h ~ Cs/U where 
17 is the angular velocity and Cs « ^fcT/p is the sound speed 
with /r as the mean particle mass and set H/h = 2. This 
value of H/h is consistent with the values used in D’Alessio 
et al. (2005). We use a single height for the disc wall by cal¬ 
culating h with the average temperature and for our fidu¬ 
cial model, H « 5.2 Rq. With this value for H, the width 
of region in the primary eclipse shadow is about the size 
the disc wall height and the width for the secondary eclipse 
is smaller than the scale-height so straight edges is a poor 
approximation. Since the shadows of the stars on the disc 
don’t strongly affect the light curve, this neglect should not 
strongly affect our conclusions. 

For ease of comparing the light curves, the flux for each 
light curve is converted into a magnitude and then the mag¬ 
nitude of each light curve is scaled so that the magnitude 
at <j) = 0 for every light curve is exactly equal to one. Since 
the magnitude is scaled, the distance to the system, d, is 
arbitrary. 

For the fiducial model. Figure 4 shows the azimuthal 
temperature profile on the left panel and the light curves 
at various wavelengths in the middle panel. The light curve 
has a smooth non-sinusoidal shape. The luminosity ratio of 
the primary to the secondary is about 3 so the peak in the 
temperature from the secondary is nearly the same size as 
the primary peak. The secondary temperature peak is large 
enough to cause a second peak in the light curve that is 
about one quarter the size of the primary peak, depending 
on the wavelength. The light curve reaches a maximum when 
the primary is closest to the centre of the observable area 
of the inner disc wall at (j) = 7r/2. The light curve reaches a 
minimum after (f> = n and again before (j) = 2tt since there is 
a smaller peak at <j) = 3 t:/2 when the secondary is close to 
the midpoint of the observable side of the disc. Despite our 
care in computing shadows in the temperature profile they 
do not significantly affect the light curves. 

The size of the light curve peak depends on the wave¬ 
length observed. For the fiducial model the size of the mag¬ 
nitude difference. Am = nimax — mmin, in M’ band (4.68 
^m) is about 0.08 mag. The peak shrinks quickly with de¬ 
creasing wavelength from 4 fj,m even though the peak wave¬ 
length for a 900 K blackbody is ~ 3.2 /rm. At these shorter 
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Figure 4. Fiducial binary and circumbinary disc model. The left panel shows the azimuthal temperature profile of the inner disc wall. 
The centre panel gives the light curves at five different wavelengths. Standard IR filter wavelengths were chosen: long dashed magenta 
for H (1.6 pm), dashed red for K (2.2 pm), solid yellow for L’ (3.77 pm), dash-dot green for M’ (4.68 pm), long dash-dot blue for N (10.5 
pm). The right panel is the peak to peak magnitude difference at different wavelengths. Parameters for the model are listed in Table 1. 


wavelengths, the light contribution from the stars is much 
larger so the variability of the disc is reduced. For wave¬ 
lengths longer than 4 pm, the light curve peak slowly de¬ 
creases with the disc wall flux. The right panel of Figure 
4 shows the wavelength dependence of the magnitude dif¬ 
ference. The wavelength for which the temperature varia¬ 
tion effect is the largest depends on the temperature of the 
disc. For the fiducial model, the light curve modulations are 
strongest in the M’ band. Hotter discs emit more strongly 
at shorter wavelength but there is also a stronger contribu¬ 
tion from the central stars at those wavelengths so the peak 
remains smaller in K than M’ except for high disc tempera¬ 
tures. The light curve variations also decrease more quickly 
with longer wavelength for hotter discs. For cooler discs, the 
largest light curve amplitude occurs at longer wavelengths. 
The temperature of the central stars also affects the wave¬ 
length dependence. For hotter stars, the light shorter wave¬ 
lengths have larger contributions from the stars and the size 
of the magnitude differences decreases more quickly. Cooler 
stars contribute less light so the variations from the disc are 
more prominent. For most low-mass binaries, the magnitude 
differences are largest between 3 pm and 10pm. 

We varied the inclination of the system from fiducial 
model value of the 75° from 40° to 89°. The inclination ef¬ 
fects the size of the light curve variations but not the shape 
significantly, see Figure 5. Some small changes in the shape 
of the light curve occur from the edge of the disc obscuring 
itself at high inclinations but the most significant change is 
in the size of the magnitude difference. The largest magni¬ 
tude difference of ~0.1 mag for the fiducial model occurs 
at about 84°. At higher inclinations, the amplitude drops 
very quickly because of the disc occultation. The inclina¬ 
tion for the maximum amplitude depends on the geometry 
of the disc. Increasing the radius allows for higher inclina¬ 
tions without occultation while increasing the height of the 
disc requires lower inclinations. For the fiducial model, the 
binary star starts to eclipse the disc and the disc eclipse the 
binary for inclination over about 83° but those effects are 
neglected here. Those occultation would have significant ef¬ 
fects on the light curve as they cause decreases in flux larger 


than variations predicted by this model. Light curve mod¬ 
ulation from variations in the temperature profile are most 
likely to be observed at inclination slightly below the limit 
for the disc occulting the stars. However, the light curve 
modulations continue to prominent at moderately low incli¬ 
nations. There is only a 0.02 mag reduction in the magnitude 
difference when the inclination is decreased from 70° to 40°. 

The light curve variations depends significantly on sev¬ 
eral parameters: stellar luminosity ratio (L*,i/L*, 2 ), stellar 
mass ratio {Mi/M 2 ), orbital separation, and radius of the 
disc. We study each dependence individually using the fidu¬ 
cial model as a starting point. 

3.1 Stellar Luminosity 

To study the effect of the stellar luminosity on the light 
curve, we kept all parameters from the fiducial model con¬ 
stant while varying the luminosity of the primary star. The 
luminosity was adjusted by changing the temperature and 
radius of star 1 in a way that is consistent with the evolu¬ 
tionary models of Siess et al. (2000). This is equivalent to 
changing the age of star 1. The parameters of star 2 are kept 
constant. We also adjusted the height of the disc to account 
for the different average temperature of the disc wall and 
to keep the hydrostatic equilibrium condition satisfied. The 
luminosity was set to 0.4Li, 0.5Li, Li, 2Li, and 5Li where 
Li is the luminosity of star 1 in the fiducial model. The 
lowest value of luminosity used is equal to the lowest possi¬ 
ble luminosity along the evolutionary track. Decreasing the 
luminosity is equivalent to increasing the age of the PMS 
star. The lowest luminosity occurs at about ~10 Myr and 
the star is half the fiducial value at about 5 Myr. Decreasing 
the age to ~0.8 Myr and ~0.3 Myr increases the luminosity 
twice and five times the fiducial value. The primary is star 1 
for every luminosity except for 0.4Li where star 1 is slightly 
less luminous than star 2. 

The luminosity ratio, L 1 /L 2 , changes the depths of the 
dips in temperature due to stellar eclipses and changes the 
shape of the SED of the stars. The amplitude of variability 
in the azimuthal disc temperature increases as the stellar 
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Figure 5. Fiducial binary and circumbinary disc model but at different inclinations. The azimuthal temperature is shown in the left 
panel and is the same as in figure 4. The centre panel shows the M-band light curve but with i = 89° in long dashed magenta, 85° in 
dashed red, 70° in solid yellow, 55° in dash-dot green, and 40° in long dash-dot blue. The right panel shows the M-band peak-to peak 
magnitude difference at different inclinations. The magnitude difference is largest at about 84°. 




Figure 6. Fiducial model but for primary stars of varying luminosities consistent with a 1 Mq PMS star. Secondary star parameters 
kept constant. The left panels show the azimuthal temperature profiles and the centre panels are the associated M’ band light curves 
for each temperature profile. Line styles and colours mark which temperature profile goes with which light curve.The lines are longed 
dashed magenta for 0.4Li, dashed red for 0.5Li, solid yellow Li, dash-dot green for 2Li, and long dash-dot blue for 5Li with Li that 
of the fiducial model. The right panels show the magnitude differences in the light curves at three wavelengths, 2.2 (K), 4.68 (M’) and 
10.5 pm (N). 


luminosity ratio departs from the fiducial model value (see 
left panel of Figure 6). The temperature decreases with lumi¬ 
nosity but the variability increases. The size of the variation 
reaches a minimum in both the temperature profile and the 
light curve when the primary to secondary luminosity ratio 
is equal to the mass ratio, Mi/M 2=2 which is slightly below 
the fiducial model luminosity ratio of 2.795. At this lumi¬ 
nosity ratio, the temperature peaks are equal which leads 
to a sinusoidal light curve with a period twice that of the 
orbital period. The effect is very similar to when there are 
two identical stars. Increasing the temperature variations 
by either decreasing or increasing the luminosity ratio from 
2 causes an increase in the magnitude difference and the 
shape of the light curve becomes more sinusoidal as shown in 
the centre panel of figure 6. The temperature dips from the 
shadow at 7 = TT increases with temperature which causes a 
slightly larger dip in the light curve. The dependency of the 


magnitude difference on wavelength is similar to the fiducial 
model except at the largest luminosity ratio. In the M’, and 
N bands, the magnitude differences are are similar and the 
difference in K band is much lower except for larger lumi¬ 
nosities. At 5Li, the magnitude difference is no longer the 
smallest in K but is the smallest instead at a much longer 
wavelength in N band. For more luminous stars, the disc wall 
temperature is much hotter so that the largest magnitude 
differences occur at shorter wavelengths in the near-IR. De¬ 
spite the unrealistic temperature adjustments to the model, 
we expected the behaviour of the light curves to be very 
similar if more the luminosity was adjusted in a way more 
consistent with stellar models. 
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Figure 7. Fiducial model but with star 2 varying mass and varying luminosity consistent with a 2 Myr PMS star using evolution models 
from Siess et al. (2000). Star 1 parameters are kept constant. The left and the centre panels are the temperature profiles and the light 
curves, respectively. Matching line styles and colours mark temperature curves with its associated light curve. The companion masses 
are 0.1 (long dashed magenta), 0.6 (dashed red), 1.0 (solid yellow), 2.0 (dash-dot green), and 3.0 Mq (long dash-dot blue). The right 
panels are the magnitude differences as a function of star 2 mass in K (red), M’ (yellow) and N band (blue). In the top panels, star 1 
mass is 1 Mq. In he bottom panels, star 1 mass is 0.5 Mq. Star 1 is not always the primary so for a larger mass star 2, the peak of the 
light curve is at </>=37r/2 instead of 0=7r/2. 


3.2 Stellar Mass 

In this section, we vary the mass of one of the stars while 
keeping the stellar parameters consistent with Siess et al.’s 
(2000) stellar evolution models of 2 Myr old PMS star. First 
we keep one star the same at 1 Mq while varying the com¬ 
panion over a range of masses (see top panels of Figure 7) 
and then keep one star at 0.5 Mq while varying the compan¬ 
ion over the same mass range (see bottom panels of same fig¬ 
ure). The companion masses range from 0.1 to 3.0 Mq. The 
orbital separation remains the same as the fiducial model at 
30 Rq for each pair and the inclination at 75°. 

The stellar mass ratio sets the semi-major axes of the 
central stars. As the mass of the star increases, it moves 
closer to the centre of mass and with a smaller orbit, the 
size of the distance oscillations shrink compared to the dis¬ 
tance to the disc wall. However, the luminosity increases 
with mass so the temperature of the disc wall increases along 
with its variations. The light curve becomes more sinusoidal 
like the previous section. With decreasing mass of star 2, 
the orbit of star 1 shrinks. With smaller distance oscilla¬ 
tions and less stellar irradiation, the disc wall cools and the 


azimuthal variations in temperature become smaller. In the 
light curve, the decrease in temperature variations leads to 
smaller magnitude difference and a less sinusoidal curve. As 
the primary moves inward with increasing mass ratio, the 
secondary moves closer to the disc wall. When close enough 
to the wall, the secondary becomes the dominant source 
of temperature variation and the peak of the light curves 
changes from n/2 to 37r/2. In the bottom centre panel of 
figure 7, this occurs between the red and yellow lines which 
are for star 2 masses of 0.1 and 0.6 Mq and star 1 mass of 
0.5 Mq. This does not occur for the 1.0 Mq star. 

The largest difference occurs for companions with 
masses larger than 2 Mq for both masses but the magni¬ 
tude differences are bigger for the 1 Mq star. There is a 
local minimum in the magnitude difference when the stars 
are identical that is more pronounced for the 1 Mq star, see 
the left panels of Figure 7. In the upper left panel, there 
is a small peak in magnitude difference at about 0.5 Mq 
companion mass for the 1 Mq star and then the difference 
steadily decreases. For the lower masses, the companion has 
too low of a luminosity to induce strong variations in the 
temperature despite being closer to the disc and the pri- 
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Figure 8. Fiducial model orbital parameters with star 2 varying mass and varying luminosity consistent with brown dwarfs at age 1 
Myr using models from Baraffe et al. (2003) and star 1 parameters are kept constant as a 1 Myr PMS star. In the top panels, star 1 
mass is 1 Mq. In he bottom panels, star 1 mass is 0.5 Mq. The left and the centre panels are the temperature profiles and the light 
curves, respectively. Matching line styles and colours mark temperature curves with its associated light curve. The companion masses 
are 0.01 (long dashed magenta), 0.02 (dashed red), 0.04 (solid yellow), 0.06 (dash-dot green), and 0.08 Mq (long dash-dot blue). The 
right panels are the magnitude differences as a function of star 2 mass in K (red), M’ (yellow) and N band (blue). 


mary’s orbit becomes too small to significantly vary the dis¬ 
tance between it and the wall. The 0.5 Mq does not have a 
sharp decrease in magnitude difference for an identical PMS 
companion and stays around 0.8 mag in M’ band within the 
range 0.3 to 1.0 Mq. 

To explore masses lower than 0.1 Mq, we change the age 
of the binary star to 1 Myr and use brown dwarf evolution¬ 
ary models from Baraffe et al. (2003) for the companions. 
The primary is again a 1 Mq PMS star and then a 0.5 Mq 
PMS star. The mass of the secondary is varied from 0.01 to 
0.1 Mq which covers the brown dwarf companion regime. 
The temperature profiles, light curves and magnitude dif¬ 
ferences are shown in the left, centre and right panels in 
Figure 8, respectively. For both primary masses, the mag¬ 
nitude difference decreases with decreasing secondary mass 
but the decline is quicker for the 0.5 Mq star. As the system 
ages, the magnitude difference will remain about the same 
as the brown dwarf cools and retracts since the PMS star is 
the source of most of the light curve variations. The dip in 
the temperature profile and the light curve decrease as the 
eclipse shadow region decreases but most of the light curve 
variation is caused by the primary unevenly irradiating the 
disk so the effect of decreasing eclipse shadow is small. Low 


mass companions cause more variability in the light curve 
for low mass primaries. 

3.3 Binary Separation 

To study the dependency on the binary separation, we vary 
the separation while keeping most of the fiducial parame¬ 
ters constant. Since the disc wall radius is defined in re¬ 
lation to the orbital separation, varying the separation re¬ 
quires adjusting the radius of inner disc wall. The height 
also increases with the disc radius even though the disc wall 
temperature decreases. The stellar properties are kept to the 
values in the fiducial model. In Figure 9, we vary the orbital 
separation from 10 Rq to 120 Rq. As the orbital separa¬ 
tion increases, the temperature of the disc wall decreases 
and the variations in temperature decrease proportionally. 
The shape of the light curve remains the same with varying 
orbital separation. Only the peaks change in size, increas¬ 
ing with decreasing temperature. The magnitude difference 
is very colour sensitive. The amplitude of the light curve 
variations increases with increasing separation in M’ and N 
band but in K, the amplitude peaks about 0.04 mag around 
20Rq . This occurs because blackbodies peak in the K band 
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Figure 9. Fiducial model but with changing orbital separation. The left and centre panels are the temperature profiles and the light 
curves, respectively, with matching line styles matching profiles with light curves. The separations used are 10 Rq (0.046 AU), 30 Rq 
(0.14 AU), 60 Rq (0.28 AU), 90 Rq (0.42 AU), and 120 Rq (0.56 AU) in long dashed magenta, dashed red, solid yellow, dash-dot green, 
and long dash-dot blue. The fiducial model is the dashed red line. The right panel is the magnitude difference versus orbital separation 
in K (blue), M’ (yellow), and N band (red). 


at about 1300 K and the maximum variation occurs near 
the peak wavelength of the blackbody. For the M’ band, 
the temperature is 600 K so the magnitude difference starts 
to decrease around 120 Rq. For N, the temperature for the 
peak wavelength to be that band is about 300 K so the mag¬ 
nitude difference does not reach its maximum in that band 
in the range we studied. Binaries with longer periods have 
larger light curve modulation in the mid-IR. 


3.4 Disc Radius 

The last parameter we explore is the circumbinary disc ra¬ 
dius. In the fiducial model, the disc radius is set to the min¬ 
imum dynamically stable orbit, Ro = l.SOaep = 54 Rq. The 
disc wall can be much farther from this point because of 
processes such as dust sublimation or photoevaporation. A 
strong photoevaporation wind can slightly increase the disc 
wall radius (Alexander 2012). We increase the radius by a 
factor of 1.5, 2, and 2.5. Generally, the temperature profiles 
flatten with larger disc radii with the dips in temperature 
from the stellar eclipse shadows decreasing slowest (see left 
panel of Figure 10). Unlike the orbital separation case, the 
light curve shape changes significantly with increasing disc 
radius. Both peaks in the light curve shrink and the smaller 
peak becomes a dip (see centre panel Figure 10) and resem¬ 
bles dips caused by extinction. The peaks in the temperature 
curve flatten which leaves only the temperature dips from 
the stellar eclipses as the source of variation. At all wave¬ 
lengths studied, the magnitude differences decrease with in¬ 
creasing wall radius with the longer wavelengths decreasing 
more slowly. The light curve variations in K reach 0 at about 
Rd = 140Rq. The differences between light curve variations 
between M’ and N are small (<2) in the radius range. 

Dust sublimation occurs at approximately 1500 K. For 
the fiducial model, the disc is already well below this tem¬ 
perature but for a closer binary, the temperature can be 
above this limit if the disc wall is at the minimum stable ra¬ 
dius, see 10 Rq separation profile in Figure 9. For binaries 
closer together than our fiducial model, there is an initial 


Table 2. Coku Tau/4 Parameters 


Ref. 

Ml 

0.6 Mq 


1 

M2 

0.5 Mq 


1 

Ri 

1.261 Rq 


1 

R2 

1.178 Rq 


1 

Teff.l 

4027 K 


1 

Teff,2 

3894 K 


1 

^sep 

1598 Rq ? 

e7.43 AU 

2 

Rd 

2877 Rq ? 

el3 AU 

2 

i 

60 deg 


2 

H/Rd 

0.11 




Parameters for model of Coku Tau/4 system 
References; 1- Siess et al. (2000) ; 2- Nagel et al. (2010) 

increase in the light curve amplitude as the temperature of 
the disc decreases with increasing radius before the ampli¬ 
tude decreases like in Figure 10. This initial increase is from 
the peak blackbody wavelength increasing to the IR bands. 
Generally, discs with wider clearings have much smaller vari¬ 
ations. 


3.5 Application to CoKu Tau/4 

We apply our model to the known binary system CoKu 
Tau/4 that Nagel et al. (2010) predicted would have light 
curve variations from the temperature oscillations. The 
Goku Tau/4 system is binary star comprising similar mass 
PMS companions and a circumbinary disc. To apply our 
model to this system, we use the same stellar parameters 
used by Nagel et al. (2010) for their SED model. The age 
of the stars is set to 4 Myr and the properties of the stars 
are taken from evolution models by Siess et al. (2000). We 
use Osep = 1598 Rq = 7.43 AU, an inclination of 60° and 
Rd = l.Sosep which are the same values Nagel et al. (2010) 
used for their fiducial model. We use a slightly different value 
for the height of the disc wall. Our wall height is set by hy¬ 
drostatic equilibrium like our fiducial model which is slightly 
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Figure 10. Fiducial model but with changing radius of the disc wall. The left and centre panels are the temperature profiles and light 
curves with matching line styles and colours. The disc radii used are 54 Rq {Rd, fid = 0-25 AU), 81 Rq (l.SiJ^i au = 0.38 AU), 108 Rq 
(2i?£i Au = 0.5 AU), 135 Rq (2.5i?i3 aa = 0.63 AU) in long dashed magenta, dashed red, solid yellow, and dash-dot green. The right 
panel is the magnitude difference in K (blue), M’ (yellow), and N band (red). 
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Figure 11. Light curve model for the CoKu Tau/4 at different wavelengths. Left panel shows the azimuthal temperature profile. The 
middle panel is the light curves at 6.0 fira (long dashed magenta), 10.5 fim (dashed red), 15.0 fira (solid yellow), 20.0 fim (dash-dot 
green), and 24.0 fim (long dash-dot blue). The right panel is the magnitude difference vs wavelength. 


smaller than the height used by Nagel et al. (2010). All the 
model parameters are list in table 2. We computed the light 
curve of the system within the wavelength range of 5 fim to 
30 fim (see Figure 11). 

The light curve has two peaks of similar height since the 
stars have very similar luminosities and masses. We find no 
variation in the light curved for wavelengths less than 9 fim 
but predict ~10% variation at 10 fim and about a 20% am¬ 
plitude variation for wavelengths larger than 15 fim. Nagel 
et al. (2010) predicted ~10% variation in the mid-IR SED of 
the CoKu Tau/4 system over half it’s orbital period which is 
in good agreement with our model. Despite the strong 10/rm 
silicate feature in the spectra of their model, our model’s pre¬ 
diction of size of the variations is within about a factor of 
a couple. The largest difference between the models is our 
model predicts a peak magnitude difference of 0.28 mag at 
15 fim that does not appear in their spectra. Their model 
predicts an approximately 0.1 mag difference over the en¬ 
tire range studied with little dependence on the wavelength 
(see their Figure 11). The dust feature extends to about 15 


fim and likely diminishes the temperature oscillation effects 
near our predicted peak. We expect the inclusion of dust 
features to change the wavelength dependence of the ampli¬ 
tude of the light curve variations but should not significantly 
change the shape. 

3.6 Other Sources 

There are multiple sources of IR variability that likely occur 
simultaneously Cody et al. (2014). Variability in the near 
and mid-IR is most commonly attributed to star spots, ac¬ 
cretion, and disc extinction (e.g. Wolk et al. 2013; Parks 
et al. 2014). Accretion variations should not be periodic and 
the time-scale depends on the mechanism causing the vari¬ 
ation (e.g. Carpenter et al. 2001). Disc extinction can be 
periodic or occur on long time-scales. Spots are periodic on 
the timescale of stellar rotation. 

We listed periodic sources of IR variability in Table 3 
with estimates of the magnitude variance. The cool and hot 
spots were calculated with a simple spot model Carpenter 
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Table 3. IR Variability Sources 


Source AM’ 


Disc Wall Temperature Modulations 0.08 
Cool Spotsi 0.2“/0.07*' 

Hot Spotsi 0.4“/0.l'’ 

Disc Occultation'’ 0.1 

Stellar Eclipse^ 0.5 


Notes: a- Spot covers 30% of star; b- Spot covers 10%; c- Refers 
to primary eclipse 

References: 1- Carpenter et al. (2001) 

et al. (2001) using a spot that is half and twice the temper¬ 
ature of the fiducial primary, respectively, covering 30% of 
the star and covering 10% of the star. Cool spots variations 
are about the same size as the variations from the temper¬ 
ature modulations but the spot light curve has a sinusoidal 
shape instead of double peaks of this model. Hot spot are 
a factor of a few larger and the light curve, if the spots are 
long lived, is the same as the cool spots. The stellar eclipse 
calculated is the primary eclipse of a binary set the fiducial 
model stellar parameters and no circumbinary disc. The disc 
occultation is also calculated with the fiducial model stellar 
parameters and orbital separation. Table 3 list the size of 
the primary stellar eclipse and the size of the dip if the disc 
completely occults the primary. The disc wall temperature 
modulations are slightly smaller than the other sources of 
periodic variations but the light curve shape is unique. 


4 CONCLUSION 

Temperature oscillations along the inner wall of a circumbi¬ 
nary disc arise from the phase dependence of the distance 
between each star in a binary and positions along the disc 
wall. These temperature variations can give significant pe¬ 
riodic modulation in a photometric infrared light curve. We 
compute a magnitude difference greater than 5 per cent for 
wavelengths longer than ~ 3 pm in a light curve from a pre- 
main sequence solar mass binary with mass ratio 2, period of 
15 days and inclination of 75° embedded in a circumbinary 
disc with a cylindrical inner edge. The maximum magnitude 
difference was about 10 per cent at ~ 4 pm. We find that 
the light curve exhibits a non-sinusoidal shape consisting of 
two peaks of different heights. One star can block the light of 
the other, shadowing the disc wall, however, these shadows 
do not cause significant variations in the integrated light 
curve. The size and shape of the light curve modulations 
are sensitive to the primary to secondary luminosity ratio, 
the mass ratio, binary orbital separation and the radius of 
the circumbinary disc wall. Larger luminosity or mass ratios 
can increase the size of the estimated photometric variations 
by a factor of a few and shape of the light curve becomes 
more sinusoidal. Increasing the orbital separation increases 
the size of the light curve variations at wavelengths longer 
than ~ 2.5 pm and does not affect the light curve shape. 
Increasing the disc radius decreases the magnitude differ¬ 
ence and changes the shape from two peaks to a small peak 
with a large dip. Light curve modulations from temperature 
oscillations are strongest for two dissimilar mass stars with 
orbital periods of the order of a couple weeks at wavelengths 
over ~ 3 pm. 


Because most young stars are born as binaries, infrared 
photometric surveys of young stars should find a class of 
periodically varying objects that are associated with emis¬ 
sion from circumbinary disc walls. The shapes of the light 
curves we have illustrated here may aid in identifying pe¬ 
riodic modulations associated with circumbinary material. 
Significant amplitude variations could be present even when 
a binary is oriented with an inclination as low as 40° to the 
line of light, so a binary need not be eclipsing to be identified 
if a circumbinary disc is present. While we have restricted 
this initial study to cylindrical shaped disc walls and cir¬ 
cular binary orbits, more complex models might make it 
possible to constrain the shape of a circumbinary disc edge. 
The photometric level of accuracy required to detect mod¬ 
ulations due to a circumbinary disc is not so exacting as to 
require a space observatory. In comparison reflection, and 
relativistic beaming require a much higher level of photo¬ 
metric accuracy. Beaming studies have discovered hundreds 
of non-eclipsing binary systems (Tal-Or et al. 2015). Hence 
infrared photometric surveys (e.g. Cody et al. 2014) may be 
quite sensitive to previously undetected circumbinary discs. 

Our model does not include any radiative transfer and 
emission from disc beyond the circumbinary disc wall was 
ignored. Depending on how the fluctuations in temperature 
along the disc wall propagate outward, radiative transfer 
could strengthen or diminish the amplitude of light curve 
variations. Demidova et al. (2013) found that the radius to 
which illumination asymmetries created by shadows of small 
mass companion can be observed depends on the viscosity 
of the disc. Asymmetries are more quickly reduced in low 
viscosity or ’’cold” discs but such asymmetries only extend 
out a small fraction of the total disc. However, if the disc 
emissions are dominated by the wall, then our model gives 
a good approximation for the shape of the light curves and 
the size-scale of flux variations. 

Thermal time-scales were also ignored in our model. 
If the time for thermal equilibrium is much shorter than 
the orbital time-scale, then neglecting heating and cooling 
time-scale is a good approximation. However if the time- 
scale is longer or comparable to the angular velocity of the 
disc wall, this temperature oscillations would smear out and 
there would be no light curve variations. If there is a slight 
lag in the thermal equilibrium, this could create an offset 
between the light curve peaks and the stellar eclipses. 

The thin cylindrical shape for the inner edge of the cir¬ 
cumbinary disc is a unrealistic simplification. An eccentric 
disc wall would impact the temperature along the edge and 
the observed flux modulations. The strength of the effects 
from an eccentric disc wall should depend on the orientation 
of the ellipse and its eccentricity. The more distant regions 
of the wall would be much cooler and vary less from star to 
disc distance oscillations. Eccentric binary orbits would also 
affect the flux modulations and depend on the orientation 
of the orbits since some regions of the disc wall would vary 
more in temperature from the stars’ movements. 
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